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ABSTRACT 

The role of radiative cooling during the evolution of a bow shock was studied in laboratory- 
astrophysics experiments that are scalable to bow shocks present in jets from young stellar objects. 

The laboratory bow shock is formed during the collision of two counter-streaming, supersonic plasma 
jets produced by an opposing pair of radial foil Z-pinches driven by the current pulse from the MAG¬ 
PIE pulsed-power generator. The jets have different flow velocities in the laboratory frame and the 
experiments are driven over many times the characteristic cooling time-scale. The initially smooth 
bow shock rapidly develops small-scale non-uniformities over temporal and spatial scales that are 
consistent with a thermal instability triggered by strong radiative cooling in the shock. The growth 
of these perturbations eventually results in a global fragmentation of the bow shock front. The for¬ 
mation of a thermal instability is supported by analysis of the plasma cooling function calculated for 
the experimental conditions with the radiative packages ABAKO/RAPCAL. 

Keywords: laboratory-astrophysics, thermal instability, bow shock, counter-streaming jets, young stel¬ 
lar objects, internal shocks, radiative cooling 


1. INTRODUCTION 

One of the characteristic features of protostellar jets 
is the presence of shocks. They can be seen as large- 
scale terminal bow shocks or working surfaces (known as 
Herbig-Haro or HH objects) that form as the jet interacts 
with previous jet ejections. Smaller-scale, internal shocks 
are also present which are driven by highly-variable flow 
velocity in the jet. These are formed as the flow moves 
away from the protostar and encounters and overtakes 
slower material from previous ejections, with the process 
repeating along the jet be am. Proper motion measure¬ 
ment s of the jets in HH 34jReiDurth et al.ll2002fl and HH 
111 (|Hartigan et al.l[200111 show the jet flow can reach 
peak velocities of ^^200—300 km s“^ with a typical ve¬ 
locity variation along the jet of ~40 km s“^. 

Shocks from protostellar jets exhibit complex dy¬ 
namics in which different effects such as shear, hy¬ 
drodynamic instabilities and rad i ative cooling can be 
present simultaneously (|Hartiganl 1200311 . In this work 
we are particularly interested in the effect of radiative 
cooling, as it can drastically modify the shock mor¬ 
phology. R adiative losses are strongly dependent on 
the opacity (|Drake et al.l 120061) . and if the shock re¬ 
gion is optically thin, then radiation can escape the 
shock leading to an increase in the post-shock den¬ 
sity. As the shock cools down it can be prone to 


the g rowth of thermal instabilities (lEieldl 119651 iHunted 
119701 ) which can fragment and ultimately break up the 
shock. This effect has previously been s tudied mostly 
through numerical simulatio n s (see e.g. iBlondin et aH 


(119891); IBlondin et al.1 (I199C 

). IStone & Norman ( 

1993 

), 

de Gouveia dal Pino & Benz 

(I1993D. Frank et al.l ( 

1998 

), 


Tesileanu et alJ (120081 ) and' Asahina et al.1 (j2014D ). 


In this paper we describe laboratory experiments which 
provide a complementary approach to study the effects 
of radiative cooling on the structure of the bow shocks. 
The similarity in the key dimensionless parameters char¬ 
acterizing our experiments, such as the Mach number 
and the cooling parameter, mean that the results are 
scalable to the internal shocks observed in YSO jets. We 
observe the formation of a bow shock in the flow and 
its subsequent fragmentation, consistent with the onset 
of thermal instabilities. The bow shock in the exper¬ 
iments is formed from the interaction of two counter¬ 
streaming flows, which is equivalent to observations of 
internal shocks in YSO jets from a reference frame mov¬ 
ing with the shock. An overall approach to modelling as- 
trophysical phenomena in laboratory exp e rimen ts is pre¬ 
sented in the review bv [Remington et al.l (l2006f) . while a 
recent review of the synergy between observations, the¬ 
ory and experi ments relevan t to th e studies of YSO jets 
can be found in iFrank et al.l (j2014D . 
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2. EXPERIMENTAL SETUP 

The experimental setup is represented schematically in 
Fig. [TJ Two counter-streaming, supersonic plasma out¬ 
flows are produced using two co-axial and oppositely fac- 
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Fig. 1. — Schematic experimental configuration to study the for¬ 
mation of a bow shock from the interaction between two counter¬ 
streaming jets with different relative axial velocities, represented as 
opposite vertical arrows on axis. The schematic depicts a side-on 
(radial), cut view of the system which has azimuthal symmetry. 
The dashed (red) arrows represent the path of the current that 
drives the two plasma flows. The (blue) arrows pointing into and 
out of the page correspond to the azimuthal magnetic field gener¬ 
ated by the current, which provides the driving force for the two 
outflows. The jets are surrounded into lower density plasma (yel¬ 
low regions) which moves with the same axial velocity as the jets. 
Smaller arrows in these regions represent the plasma flow direction. 
The images depict the two counter-streaming outflows (a) before 
their collision, (b) after they collide, triggering the formation of a 
bow shock moving towards the bottom foil. 


ing radial foil Z-pinches (iSuzuki-Vidal et ^l2009t) . Each 
radial foil is a metallic disc (40 mm diameter, 14 /rm 
thick aluminum), subjected to a fast-risin g electrical cur¬ 
rent pulse from the MAGPIE generator ( Mitchell et alJ 
llQQOf l. which for the present experiments used a peak 
current of MA in ^330 ns. The current is driven to 
the foils through 6.35 mm diameter stainless steel tubes 
touching each foil at its centre, with the same current 
going through both of the foils via vertical posts (shown 
schematically in Fig. [1^.). The distance between the two 
foils surfaces was ~30 mm. 

The plasma is produced by continuous ablation of the 
surfaces of the foils as they are heated by the current, 
and the ablated plasma is accelerated by the axial pres¬ 
sure gradient produced by the current-induced azimuthal 
magnetic held as it diffuses through the foils. Each of the 
foils produces a supersonic plasma outflow which consists 
of a dense central jet surrounded by lower-density ambi¬ 
ent plasma. Both components propagate with the same 
axial velocity of ^50—100 km s“^. Details on the forma¬ 
tion of t he outflows by a single rad i al foil Z-pinch can b e 
found in ISuzuki- Vidaj et 311(1200911 . ICiardi et al.l (1200911 . 

IGourdain et al.l(|201Clll and lSuzuki-Vidal et al.l(|2012ll ~In 

the experiments presented here, the interaction of the jets 
was diagnosed from the side-on (radial) direction using 
an optical framing camera (Invisible Vision UHSi 12/24) 
that imaged the optical self-emission from the plasma. 
This camera is capable of taking up to 12 images per ex¬ 
periment, with 5 ns exposure and 30 ns inter-frame sep¬ 
aration. We also used simultaneous optical laser shad- 
owgraphy and interferometry (A=532 nm, pulse dura¬ 
tion 0.3 ns). The latter diagnostic was used to measure 
the electron density distribution of the different plasma 
features present, i.e. jets, surrounding plasma and bow 


shock (|Swadling et al.l[2014f) . 

2.1. Scaling 

Jets from young stellar objects are generally well de¬ 
scribed by ideal magneto-hydrodynamics (MHD), and 
the experiments presented here are designed to pro¬ 
duce flows in a similar regime. It is evident that these 
evolve over hugely different length and time-scales, and 
have different physical characteristics in terms of den¬ 
sity, temperature and chemical composition. Neverthe¬ 
less, the invariance properties of the ideal MHD equa¬ 
tions provide a framework that allows meaningful scaling 
of t he jet dynami c s over many orde rs of magnitude (see 
e.g. iRvutov et al.l (|1999L 120001 l2001h l. provided that cer¬ 
tain constraints between the flow variables are satisfied. 
It is important to note that such invariance is also appli¬ 
cable to flows with shocks and to r adiative flows und er a 
more restrictive set of constraints (|Falize et al.ll20li[l . 

In order for a fluid description to be applicable, we re¬ 
quire the flows to have a localisation parameter J <|; 1. 
This is equivalent to an ion mean free path much less 
than the characteristic spatial scale of the system. U sing 
the parameters from Table 1 in iHartigan et al.1 (I2009fl we 
estimate jets from YSOs have 6 ^ 10“®. Applicability 
of an MHD description also requires that the transport 
of momentum, magnetic held and thermal energy occurs 
predominantly through advection with the how, i.e. neg¬ 
ligible dissipation through viscosity (Reynolds number 
Re ;§>!), magnetic diffusivity (magnetic Reynolds num¬ 
ber Rcm 1) and heat conduction (Peclet number 
Pe :$> 1) respectively. Due to their large spatial scales, 
YSO jets are characterised by Re, Rcm, Pe > 10^. Our 
experiments are character ised by J ^ 10~^, Re ^ 10®, 
Pe ~ 10® and Rcm^ ~ 10® ijSuzuki-Vidal et £il]l2012t) and 
thus we expect a similar overall physical behaviour. 

Additionally, the interaction of radiative jets with the 
interstellar medium can be broadly classihed by three di¬ 
mensionless parameters (|Blondin et al.l[T990ll : the Mach 
number M (the ratio of how speed to sound speed), 
the density contrast r] (the ratio of density between the 
how and the ambient where it propagates), and the cool¬ 
ing parameter Xcooi (the ratio of the cooling time TcooI 
to the characteristic hydrodynamical time of the how 
Thy dr o), which quantihes the effect of radiative losses in 
the plasma. The specihc values of M, rj and Xcooi es¬ 
sentially determine the overall morphology of the how. 
Internal bow shocks in YSO jets are characterised by a 
relative velocity of ^40 km s~ ^ and M ^ 10, and w e can 
expect ry 1 and Xcooi ^ 1 (|Hartigan et al.l[2009f) . As 
will be discussed later in the paper, the values of these 
three parameters in our experiments are close to those in 
YSO jets. 

Overall, the similarity of the dimensionless parame¬ 
ters allows applying the Eulerian s caling relations de - 
scribed in detail in iRvutov et al.l (119991 . 120001 I2001D . 
Flows with identical Mach numbers will evolve with iden¬ 
tical morphology, but on different hydrodynamic tempo¬ 
ral {Thydro) and spatial {rjet) scales, related via the cor¬ 
responding how velocities (Vfiow) as Thydro = rjet/Vfiow 
Taking the jet radius as a characteristic spatial scale, 
young stellar jets typically have ryso 50 AU ^ 
10^® mm, and experiments rexp 1.5 mm), how veloci¬ 
ties of Vyso = 40 km s“® (i.e. typical velocity variability 
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Fig. 2.— Counter-streaming jet interaction results from optical self-emission of the plasma obtained from the same experiment. The 
arrows in the first three frames indicate the position of the tip of both jets (visibility dependant on image contrast levels), with their 
collision highlighted at 370 ns. The last six frames are focused on the bow shock region, which is seen to fragment most evidently in the 
last three times. 


in YSO jets) and for the experiments here Kxp = 140 
km s“^ (i.e. the relative velocity of a single jet in our ex¬ 
periments, see explanation in the next section). Thus the 
characteristic temporal scales are tyso =1.9x10® s (^6 
years) and = 10 ns respectively. The total time in¬ 
terval over which the evolution of the flows is followed in 
the experiments of ^300 ns (« 30 Texp) thus corresponds 
to ~180 years of evolution for the astrophysical counter¬ 
part, which exceeds the typical time scale for multi-epoch 
observations of YSOs with the Hubb l e Spa ce Telescope 
of ^10 years (see e.g. iHartigan et all (1201111 b 

3. EXPERIMENTAL RESULTS AND DISCUSSION 

The overall evolution of the interaction of the two 
counter-streaming plasma jets can be seen in Fig. ^ 
which shows side-on optical emission images obtained in 
the same experiment at different times after the start 
of the current pulse driving the jets. Images in the top 
row of Fig. [2] show the entire region between the two 
foils and correspond to the early times of the interac¬ 
tion (310—430 ns). The images show the formation of 
two well-collimated jets on the axis of the system prop¬ 
agating towards each other. Although the top and bot¬ 


tom foils are nominally the same and are driven by the 
same current, an asymmetry between the top and bot¬ 
tom flows is evident from the emission images. The top 
jet is formed first compared to the bottom jet, which is 
seen by looking at the positions of the tips of the jets in¬ 
dicated by horizontal arrows on the first 3 frames. This 
asymmetry is caused by opposite polarity current drive 
for each foil (radially outwards on the bottom foil com¬ 
pared to radially inwards on the top foil) and is repro¬ 
ducible from experiment-to-experiment. The asymmetry 
affects the details of the plasm a ablation and initial acce l- 
eration of the flow at the foil (iGourdain fc Sevl^l2013[l . 
This results in different flow parameters (e.g. ram pres¬ 
sures) which, as the jet collide, lead to the formation 
of a bow shock. This is first evidenced at 370 ns as a 
highly-emitting region aligned with the head-on collision 
between the two jet tips. The bow shock is fully formed 
at 400 ns and is seen to move downwards, i.e. towards 
the bottom foil. At this time the post-shock region is 
seen as a highly-compressed, highly-collimated column, 
while the bow shock begins to develop small-scale struc¬ 
tures from here onwards. The long-term evolution of the 
bow shock and of these structures is highlighted in the 
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subsequent images, corresponding to times 460—640 ns. 

The first collision between the two jets seen in the im¬ 
age at 370 ns occurs ^2 mm below the mid-plane, which 
suggests that the two jets have slightly different veloci¬ 
ties. This is confirmed by following the positions of the 
visible tips of the jets in time before their collision, re¬ 
sulting in tip velocities of Vtop jet 80 ± 10 km s ^ and 
Vbottom jet 60 ± 10 km s“^. The higher velocity in 
the top jet, combined with a larger density inferred from 
the stronger self-emission and from laser interferometry, 
imply an imbalanced ram pressure {Pram = pV'^) in the 
collision, which explains the shape and orientation of the 
bow shock and its downward propagation. The velocity 
of the leading edge of the bow shock was measured as 
Vbow ^ 40 ± 10 km s“^. 

The velocity of the bow shock measured in the 
experiments can be compared with the velocity expected 
from a one -dimensional, m o mentu m flux con se rvatio n 
argument (iNorman et afl (1198311: iHartiganl (1198911: 


iBlondin et ah (1199011: Ide Gouveia dal Pino fc Benz 


1994fl : Ide Gouveia Dal Pinol (1200511 : iNicolai et al. 
(I2008f)l. This is convenient to do in the reference frame 


of a single jet driving a shock through a stationary 
medium ahead of it. The velocity of a bow shock Vbow 
is related to the velocity of the jet Vjet and the density 
contrast p between the jet and the pre-shock external 
medium by Vbow ~ Vjet{^+V~^^^)~^■ The measured 
experimental jet tip velocities of Vjet top '^SO km s“^ 
and Vjet bottom ~60 km s“^ for the top and bottom 
jets respectively can be translated into a single jet with 
a relative tip velocity of Vjet rei=140 km s“^ which 
in turns forms a bow shock with a relative velocity of 
Vbow rei=100 km s~^. Using these relatives velocities in 
the equation above result in a density contrast rj 6, 
i.e. the experiment is equivalent to the interaction of 
a single jet which is 6 times denser than the ambient 
medium where the shock is generated which, for this 
particular counter-streaming geometry, the ambient 
medium is the opposite (bottom) jet. We can cal¬ 
culate the internal Mach number of the jet working 
surface by taking the ratio of the jet relative velocity 
{Vjet r-ei=140 km s“^) and the typical ion acoustic speed 
in th e jet flow of Cs ~15 km s“^ (jSuzuki-Vidal et id] 
[Ml), resulting in an internal Mach number of M 10. 
Therefore the values of p and M in the experiments 
are similar to those in YSO jets (jHartigan et al.l l2009fl 
and thus the experiments should evolve with a similar 
jet/shock morphology as discussed earlier. 

In addition to the bow shock, the images in Fig. [2] also 
show the collision between the two flows off-axis, i.e. the 
plasma surrounding the jets, which leads to the forma¬ 
tion of a double-shock structure that extends at large 
radii. This feature is especially evident at the top of the 
images from ^460 ns onwards. These standing shocks re¬ 
main approximately symmetric with respect to the mid¬ 
plane between the two foils. Preliminary MHD simula¬ 
tions indicate that these features arise due to the pres¬ 
ence of toroidal magnetic held advected b y the counter¬ 
streaming flows (|Suzuki-Vidal et aI]l2014D . A more de¬ 
tailed discussion of the pile up of the magnetic fields will 
be published separately. 


3.1. Bow shock and jet working surface: evolution and 
small-scale structures 

The emission images in Fig. |2l show in detail the evo¬ 
lution of the bow shock from its formation at 400 ns 
through its rapid fragmentation. At 400 ns the surface 
of the bow shock and the post-shock region just behind it 
are smooth. The formation of small-scale structures oc¬ 
curs on timescales that are shorter than the inter-frame 
separation for this particular diagnostic. This gives an 
upper limit for their development timescale of ~30 ns, 
which is consistent with the appearance of new struc¬ 
tures between frames at later times. The characteristic 
spatial scale of the structures seen in the images ranges 
from ^200 ^m to ~2 mm, and at the final stages of the 
evolution (the last two panels at 610 ns and 640 ns) the 
bow shock has fragmented into a large number of emit¬ 
ting clumps. The smallest observed size ('^200 ^m) seen 
on the images obtained with the optical framing camera 
is comparable with the spatial resolution of this diagnos¬ 
tic, and is also comparable with the motional blurring 
due to the temporal resolution of the diagnostic (5 ns x 
40 km s“^ = 200 pm). 

Further details of the small-scale structures present in 
the bow shock and post-shock region were obtained with 
laser probing (0.3 ns pulse duration), which has signifi¬ 
cantly better spatial resolution (< 50 pm) and reduced 
motional blurring (~ 10 ^m). The interferometry chan¬ 
nel of this diagnostic provides measurements of the spa¬ 
tial distribution of the plasma electron density, and the 
shadowgraphy channel, which is sensitive to spatial gra¬ 
dients of the electron density, provides information on 
the characteristic spatial scales of the non-uniformities. 
Fig. [3] shows laser probing images obtained at 400 ns 
in the same experiment as the optical emission images 
shown in Fig.|2[ i.e. the timing of the laser probing images 
corresponds to the fourth panel in Fig.|2l Fig. [3^ shows 
the raw image from the interferometry channel. The ap¬ 
parent distortion of the interference fringes, which were 
initially horizontal and uniformly spaced, is caused by 
changes in the interference state which are induced by 
a spatially varying phase delay imparted on the probe 
beam by the pl asma. This phase del ay can be extracted 
from the image (jSwadling et al.ll2f)l^ and is proportional 
to the electron column density of the plasma, UeL (in 
cm“^). Here L is the length of plasma along the prob¬ 
ing beam, which changes as a function of position in the 
plane of the image. The sufhciently good axi al symmetry 
of the object allows applying Abel inversion ([Hutchinsonl 
l2005f) to the electron column density, resulting in the 
axi-symmetric (radial) distribution of electron density 
ne{r) (in cm“^) shown in Fig. [3)3. The highest elec¬ 
tron density near the axis, reaching rig ~10^® cm“^ (up 
to 75% uncertainty due to errors in the choice of the 
central axis and left-right asymmetries), is observed in 
the narrow, compressed region formed above the bow 
shock (the post-shock region), while the electron density 
in the flow below the bow shock (the pre-shock region) is 
He ^(3±2)xl0^® cm“^ (i.e. ~60% uncertainty). We note 
that, although these errors seem large, for measurements 
of electron density off-axis these errors decrease rapidly 
and can reach values <20%. The higher electron density 
present in the jet driving the bow shock is qualitatively 
consistent with the density contrast (heavy jet propa- 
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Fig. 3.— Optical laser probing of the bow shock at 400 ns (same experiment as Fig. [a. (a) Raw data from laser interferometry, (b) 
axi-symmetric electron density (rie) from analysis of (a), (c) Laser shadowgraphy, with the field of view shown schematically in the dashed 
inset in (a) and (b). The fringing on this image is an artefact present in this particular diagnostic. 


gating through ambient) inferred from application of the 
one-dimensional momentum flux conservation argument 
introduced in the previous section. 

The laser shadowgraphy diagnostic in Fig. [3lf shows 
several interesting features. Firstly, there are two dark 
regions in the top part of the image positioned on both 
sides from the axis, and just below them a dark verti¬ 
cal region on axis. These features are produced by large 
density gradients and their shape is consistent with com¬ 
pression of the jet to a diameter smaller than the ini¬ 
tial jet diameter due to converging flows driven by large 
post-shock pressures. Formation of such shocks was ob¬ 
served in numerical simulations of high Ma ch number jets 
propa gating through ambient media (e.g. iNorman et al.l 
(|1982fl ~). and the structure we see in the experiment is 
very similar to that observed in simulations of radia- 
tively cooled jets wit h a similar dens i ty con trast, e.g. pre¬ 
sented in Fig. 7c of iBlondin et al.l (1199011 for a density 
contrast of 77 = 8 , which is of the same order to the one 
inf erred for these exper iments of 77 = 6 . The simulations 
bv IBlondin et alJ (1199011 also show the presence of a sec¬ 
ond standing conical shock, positioned between the bow 
shock and the region of maximum compression which 
opens up towards the bow shock. The image in Fig. [3t 
indicates that this conical shock is also present in the 
experiment, though it is less pronounced indicating a 
smaller density gradient than that in the downstream 
conical shock. We also note the presence of a horizontal 
dark region connecting the dark features at the top of 
the image, which is consistent wit h the expected s hape 
of a Mach disk (see e.g. Fig, l a in iHartigaiJ (|1989ll and 
Fig. 7b in IBlondin et al.l (1199011 1. 

Finally, the image in Fig. [3]: shows the presence of 
small-scale structures already developing behind the bow 
shock at this early stage of the evolution. The pertur¬ 
bations appear to be elongated, predominantly in the 
direction normal to the jet axis, and the smallest de¬ 
tectable spatial scale is '^120—170 /im, which is a factor 
of ^3 larger than the spatial resolution of the diagnostic 
(^50 77 m in the diffraction limited case). We interpret 


the observed development of small-scale perturbations 
in the bow shock leading to its complete fragmentation 
later in time as being a result of strong radiative cool¬ 
ing in the post-shock plasma. The fragmentation of the 
bow shock and the development of dense clumps i n this 
regio n are associated with a thermal instability (iFieldl 
I1965I1 . which is discussed thoroughly in the next sec¬ 
tion. The fragmentation of shocks in protostellar jets 
due to thermal instabilities has b een previously s tudie d 
by numerical simulat i ons (see e.g. IBlondin et al.l (jl989t) : 
IBlondin et al.l (jl99r)H : iFrank et al.l (jl998l B. however the 
effect of this instability in YSO jets and shocks is still 
unclear. 


3.2. Radiative cooling and thermal instabilities in the 
bow shock 

The most striking result in the experiments is the 
rapid development of small-scale spatial structures in the 
bow shock and post-shock region. Our interpretation is 
that the fragmentation of the bow shock is related to 
a dynamic, local thermal instability that leads to the 
condensatio n of de i isity p e rturbations by radiati v e coo l¬ 
ing (see e.g. IFieldl (119651). iMathews fc Bregmaiil (I1978I1. 

I Fall fc ReesI (119851) . iBalbusI (I1986D and IBlondin fc Cioffil 
(jl989l )). The instability develops over the character¬ 
istic radiative cooling time, TcooI, which for optically 
thin plasmas, such as the ones in our experiments 
(|Espinosa et al.l [20151) . is given by the ratio of thermal 
energy density, U, to the radiated power per unit volume 
Prad = neniA{ni,Te) as Tcooi = U/neniA{ni,Te), where 
A{ni,Te) is the normalised cooling function (in erg cm^ 
s“^), and Ui is the ion density. Th is expression can also 
be written as (|Rvutov et al.lll99'^ : 


Tcooi [s] = 2.4 X 10 


(Z -b l)T[eV] 
Zni[cm-^]A{ni,Te) 


( 1 ) 


where Z is the average ionization in the plasma. 

When radiative cooling occurs on time scales that are 
long compared to the sound speed crossing time, per- 
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Fig. 4. — Calculated (a) cooling function and (b) average ionization for aluminum as a function of electron temperature for different ion 
densities relevant to our experiments (the legend shown in (a) applies to the entire figure), (c) Cooling time calculated using Eq. [T] (d) 
Thermal instability analysis applied to the cooling functions shown in (a). The regions below each horizontal line indicate the thresholds 
for the onset of isobaric (< 2), isochoric (< 1) and adiabatic (< —0.5) thermal instabilities. 


turbations over a region of size Xiso = CsTcooI (where Cg 
is the ion acoustic speed) tend to maintain a common 
pressure (isobaric). In this case nT ^constant and the 
cooling rate then scales as F = oc A{T)/T^. In the 
isobaric regime, if the cooling rate increases for a decreas¬ 
ing temperature, i.e. dV/dT < 0, then radiation losses 
will be even more efficient in removing energy from the 
plasma and further reducing its temperature. To main¬ 
tain pressure equilibrium with its surroundings the den¬ 
sity increases, thus further increasing the radiated power 
losses and potentially leading to a run-away condensation 
of the initial density pertur bations. Det a iled ana l ysis o f 
the instabil i ty wa s done by iFieldl (|1965D . iHunted ljl97O0 
and lBalbu^ (|1986f) . For isobaric modes, the condition for 
stability is given by 

d{logA{ni,T)) 
dilogT) < 

While Xiso sets an upper limit for the the length-scale 
of isobaric condensation, thermal conduction will sup¬ 
press short wavelength perturbations that are of the or¬ 
der of the so-called Field’s length {Xpieid = 27r[(7 — 
^)XthTcooiY^‘^), where Xth = n/ue is the thermal diffu- 
sivity and 7 is the ratio of heat capacities. The most 
unstable wavelength, valid in the regime of large wave- 
numbers (but still smaller than a critical wavenumber 
kcrit = ‘^'^^Fieid)y given by the geometric mean (|Fieldl 


[196^ A 

max — i^XFioidXiso') ^ SO that density and temper 
ature perturbations will be unstable in the wavelength 
range 

Xpield ^ A < Xiso (3) 

Although cooling in our experiments is different from 
that in astrophysics due to the differences in ele¬ 
ments and physical conditions of temperature and den¬ 
sities, for optically thin plasmas the relevant parame¬ 
ters for comparison are the dimensionless cooling pa¬ 
rameter Xcooi and the dependence of the cooling func¬ 
tion A(ni,T) on temperature. To estimate the effects 
of radiative cooling in the post-shock aluminum plasma, 
we use new cooling rates calculated with the compu¬ 
tatio n al packages ABAKO /RAPCAL ([Rodriguez et al.l 
I2008I iFlorido et al.l 120091) . The codes calculate the 
plasma level populations and average ionizations by solv¬ 
ing the set of rate equations of the collisional-radiative 
model implemented in ABAKO (assuming the plasma is 
optically thin and in steady-state). The model is capa¬ 
ble of accounting for coronal equilibrium, local and non¬ 
local thermodynamic equilibrium regimes. The ato mic 
data required were obtained using the FAC code (|Gul 
I 2 OO 8 D in the relativistic detailed configuration account¬ 
ing approach, with the spin-orb it split arrays formal¬ 
ism ([Bauche-Arnoult et al.l[l98^ and including configu¬ 
ration interaction within the same non-relativistic atomic 
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configurations. The databases of cooling functions and 
average ionizations obtained with ABAKO/RAPCAL 
were subsequently parametrized as a function of the 
plas ma density and tempe rature using the PARPRA 
code (jRodriguez et al.|l2ni4[l . Results from these numer¬ 
ical calculations are presented in Figs. |4^-b showing the 
variation of the cooling function A(ni,Te) and the av¬ 
erage ionization Z for an aluminum plasma as function 
of electron temperature (ranging from ^ 10^ — 10® K, 
~ 1 — 85 eV) for different ion densities typical of our 
experiments (jii = 10^^ — 10^® cm“®). Fig. |4^ shows 
the rate of cooling increases with decreasing ion density, 
varying up to two orders of magnitude. The average 
ionization in Fig. shows an increase with tempera¬ 
ture, with little variation as a function of ion density. 
With the cooling function and the average ionization we 
can calculate the cooling time TcooI (Eq. [TJ , and this is 
plotted for different aluminum ion densities in Fig. 0]:. 
The plot shows that overall the cooling time decreases 
with increasing ion density, as TcooI oc l/A(ni,Te), and 
longer cooling times can be expected around an electron 
temperature of Tg ~10® K. With the cooling time we 
can calculate the cooling parameter Xcooi = Tcooi/Thydro 
which quantifies the importance of radiative losses in the 
plasma. Taking the estimated hydrodynamical time of 
the experiments of Texp 10 ns (see Sec. EH) , we see 
that, independent of ion density TcooI > T^xp in the re¬ 
gion around Tg ^10® K, which implies Xcooi ^ 1 and thus 
radiative losses are not important. This region correlates 
with a dip in the cooling functions at this temperature. 
Besides this temperature region, however, the cooling 
time is overall of the order of the hydrodynamical time 
and thus the cooling parameter Xcooi 1- The hydrody¬ 
namical time can also be estimated as the time it takes for 
the shock to form and become fully disrupted, which from 
Fig. [2] it can be taken as Texp 550 — 370 ns = 180 ns 
resulting in Xcooi (3 — 30 ns)/180 ns < 1 and thus the 
plasma in the shock is expected to be radiatively cooled. 

Knowledge of the dependence of the cooling function 
as a function of temperature allows performing a thermal 
instability analysis. The stability condition for isobaric 
modes presented in Eq. ^ {dlogA(ni,Te)/dlogTe < 2 ) is 
plotted in Fig. |4ji using the cooling functions for differ¬ 
ent ion densities in Fig. |4ji. The threshold for isobaric 
instabilities shows that, almost independently of the ion 
density, the plasma is expected to become unstable at 
temperatures of Tg ^ (1.5 — 9) x 10^ K (~ 5 eV) and 
at Tg > 2.5 X 10 ^ K (> 20 eV). Following the analy¬ 
sis in iShchekinovI (|1978r) we also plot other instability 
thresholds, namely isochoric and adiabatic modes, which 
correspond to dlogA{ni,Te)/dlogTe < 1 and < —0.5 re¬ 
spectively. 

Furthermore, the onset of thermal instabilities in the 
post-shock plasma depends on the thermal evolution of 
the ion and electron plasma components. Because of 
the large difference between the ion and electron masses, 
there is a large temperature difference between the ion 
and electron components immediately behind the bow 
shock. In the strong shock approximation (Mach num¬ 
ber M ^ 1), the post-shock (PS) ion temperature 
is T,,p5 = 2(7 - l)AmpV,l^ rei/ksh + 1)' 10® K, 

where A is the atomic weight (A=27 for aluminum), 
nip is the proton mass, Vbow rei is the shock velocity in 


the reference frame of a stationary pre-shock medium 
(i.e. Vbow rei = 100 km s“^), and 7 = 5/3 assuming 
an ideal gas. Compared to the ions, the electrons are 
compressed adiabatically across the shock and their tem¬ 
perature increases only by a factor of ^ 2.5 to 

Te,ps ^ (2 — 3) X 10® K. Although the plasma can be 
unstable at these temperatures, the existence of a well- 
defined unstable range of wavelengths, as given in Eq. [2 
is not satisfied. The post-shock is characterized by a re¬ 
laxation region where the energy exchange between ions 
and electrons competes with radiative cooling. Immedi¬ 
ately after the shock the electrons are thermally decou¬ 
pled from the ions and their cooling time-scale (< 1 ns) is 
much shorter than the ion-electron energy equilibration 
time-scale (rg^® = jZ ~ 20ns). Thus, over time- 
scales of a nanosecond, the electron temperature rapidly 
drops while the ion temperature remains essentially con¬ 
stant. The cooling rate also decreases considerably and 
an equilibrium is immediately reached behind the shock 
where electron heating due to the energy received by the 
ions balances their radiative energy losses 

Z^«(7-l)^A(n„rg) (4) 

In this regime ss TcooI 15 — 20 ns and the elec¬ 
trons are approximately isothermal. Detailed calcula¬ 
tions of temperature equilibration show that the equilib¬ 
rium electron temperature in the post-shock relaxation 
layer is Tg ^ (1 — 1.5) x 10® K, thus the plasma is ex¬ 
pected to be thermally stable. After ^ 40 ns the electron 
and ion plasma components are fully equilibrated and 
their common temperature decreases below < 7 x 10"^ K. 
This places the plasma in conditions corresponding to a 
thermally unstable regime. The condition on the wave¬ 
lengths is satisfied and density and temperature pertur¬ 
bations with wavelengths A ^ 30 — 100 /rm grow over very 
short time-scales of order < 1 ns. 

4. CONCLUSIONS 

We presented a new experimental configuration that 
aims at studying experimentally the formation of bow 
shocks relevant to those present in young stellar jets. 
In our experiments we produce two counter-propagating, 
supersonic plasma jets from plasma ablation of two radial 
foil Z-pinches. A bow shock is driven by the head-on col¬ 
lision between the two jets and the bow shock properties 
are determined by their relative velocities and densities. 
A key result is that initially the bow shock is smooth 
but then quickly, within time scales of ~30 ns, develops 
small-scale spatial features behind the shock front. This 
time-scale is consistent with the expected cooling times 
in the experiments, which were estimated using the cool¬ 
ing functions and average ionization calculated with the 
radiative packages ABAKO/RAPCAL. This allowed per¬ 
forming a thermal instability analysis for isobaric modes 
resulting in expected temperature ranges at which the 
plasma should become thermally unstable. Detailed an¬ 
alytical calculations of the thermal evolution of the ion 
and electron components predict that, for the experi¬ 
mental shock conditions, we expect the development of 
an isobaric thermal instability in time scales of ^^40 ns, 
and with typical spatial scales of ^30—100 fim. This 
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Fig . 5.— Multi-epoch observations of HH 1 showing the formation of knots in the bow shock (adapted from Fig. 2 in IIHartigan et al.l 
1201111 '). The image shows two out of three epochs of HST images, with Ha in green, the [Su] in red, and yellow denotes emission in both 
filters. 


is in very good agreement with our experimental results 
that show an initially smooth bow shock which, within 
^^30 ns, fragments into small-scale spatial features with 
typical sizes of ^ 120 /im. 

Thermal instabilities may play a role in HH shocks 
provided the shock velocities are high enough to raise 
the temperature in the post-shock gas significantly above 
the peak of the cooling curve around 2 x 10^ K. The crit¬ 
ical shock velocity for the onset of cooling instabilities 
has been estimated t o be ~ 200 km s~^ in 1-D sim - 
ulations (|Smithl [I989I1 . though [Sutherland et al.l (|2003f) 
found that thermal instabilities in their 2-D simulations 
generated filaments and voids in the post-shock g as for 
shock velocities as low as 120 km s“^. However, (jlnnesl 
[TMl demonstrated that even a fairly weak pre-shock 
magnetic field stabilized shocks up to 175 km s“^. In 
HH jets, most shocks have velocities < 100 km s“^, and 
should not be affected by thermal instabilities. However, 
the strongest bow sh ocks, such as HH 1 and HH 2, have 
high ion ization lines (iBoehm et al.lll99^ and broad line 
prohles (iHartigan et al.lll987f) indicative of shock veloc¬ 
ities ^ 200 km s“^, well above the criteria for thermal 
instabilities. As shown in Fig. [SJ bright knots do appear 
along the bow shock in HH 1 on time-scales of decades. 
Using the scaling presented at the end of Sec. 12.11 the 
laboratory temporal scale of 30 ns would correspond to 
an astrophysical time scale of ^ 17 yr, consistent with 
the astronomical observations. The spatial scales asso¬ 
ciated with the astronomical knots also scale well with 
the non-uniformities in the experiment. For example, 
the minimum spatial scale of the non-uniformities in the 
experiments of 0.12 mm scale to 4 AU in the astro¬ 
nomical observations, where the spatial resolution is ~ 
20 AU. The larger non-uniformities present in the op¬ 
tical self-emission (e.g. the last three panels of Fig. [2l) 
have a typical size of ^ 1 mm, which scales to 30 AU in 
Fig. 5, consistent with the size of the new knots in HH 1. 
Driven from the same source on the other side of the 
outflow, HH 2 also has a high shock velocity and shows 
small knots t hat appear and merg e along the strongest 
shock fronts (|Hartigan et al.ll^Ollll . Of course, there are 
other ways to generate clumps along bow shocks such 
as Kelvin-Helmholtz instabilities or even a clumpy pre¬ 
shock density, and we cannot rule out these possibilities 
without further observational data. 

More generally, the thermal instability analysis 
could be applied to cooling curves used for nu¬ 


merical simulations of proto s tellar j ets and shock s 
(see e.g. iDalgarno fc Mc Cravl (|1972ll : iKafatosI (|1973ll : 
[Sutherland fc Dooital ( 199311 1. Although in these cases 
we expect the temperature ranges for the onset of ther¬ 
mal instabilities to differ from those expected in the ex¬ 
periments due to the different elements and abundances 
that characterize them (e.g. H in simulations compared 
to Al in the experiments), in both cases the instability 
should develop at the appropriate slope of the cooling 
curve by the condition given in Eq. [2] 

For typical interstellar medium compositions, shock ve¬ 
locities of the order of > 120 km s“^ are required to push 
the cooling curve into the regime where it declines rapidly 
enough with increasin g temperature to become prone to 
therm al instabilities (jSutherland fc Dooital 119931 : linnesi 
Il992t) . Hence, most shocks in jets do not fall into this 
regime. However, some, like the clumps observed in 
HH 2, have shock v elocities of at least 200 km s“^ 
(|Hartigan et al.l[2011ll . thus we could expect thermal in¬ 
stabilities to be related to the appearance of spatial in¬ 
homogeneities in the shock. 

The experiments presented here are important as they 
shed light on the onset and non-linear evolution of the 
thermal instability. Numerical simulations of thermally 
unstable flows can be challenging as thermal conduc¬ 
tion must be explicitly included to suppress the growth 
of small-scale perturbations. In particular, the Field’s 
lengt h has to be resolved by at le ast a few computational 
cells (|Kovama fc Inutsukall2004f) to avoid the growth of 
perturbation and fragmentation at the grid scale. Fur¬ 
thermore, the presence of magnetic fields further in¬ 
creases the complexity of numerical calculations by mak¬ 
ing thermal conduction anisotropic. In that direction, 
similar experiments to the ones presented here can be 
designed to increase sufficiently the magnetic field in the 
flow to allow studying such regime in the future. This 
could be relevant to previous theoretical studies which 
predict the sup pression of th erm al instabilities by ma g- 
netic helds (e.g. linne'sl (II992I1 and lLesaffre et al.) (1200411 b 
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